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ABSTRACT 

We present a clustering analysis of ^60,000 massive (stellar mass M* > 10 11 Mq) 
galaxies out to z = 1 drawn from 55.2 deg 2 of the UKIRT Infrared Deep Sky Sur- 
vey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II Supernova Survey. Strong 
clustering is detected for all the subsamples of massive galaxies characterized by dif- 
ferent stellar masses (M* = 10 110_11 - 5 M Q , 10 11 - 5_12 - M Q ) or rest-frame colors (blue: 
U — V < 1.0, red: U — V > 1.0). We find that more mature (more massive or redder) 
galaxies are more clustered, which implies that the more mature galaxies have started 
stellar-mass assembly earlier within the highly-biased region where the structure for- 
mation has also started earlier. By means of halo occupation distribution (HOD) 
models fitted to the observed angular correlation function, we infer the properties 
of the underlying host dark halos. We find that the estimated bias factors and host 
halo masses are systematically larger for galaxies with larger stellar masses, which is 
consistent with the general agreement that the capability of hosting massive galaxies 
depends strongly on halo mass. The estimated effective halo masses are ~ 10 14 M Q , 
which gives the stellar-mass to halo-mass ratios of ~ 0.003. The observed evolution of 
bias factors indicates rapid evolution of spatial distributions of cold dark matter rela- 
tive to those traced by the massive galaxies, while the transition of host halo masses 
might imply that the fractional mass growth rate of halos is less than those of stellar 
systems. The inferred halo masses and high fractions of central galaxies indicate that 
the massive galaxies in the current sample are possibly equivalent to central galaxies 
of galaxy clusters. 

Key words: galaxies: clusters: general - galaxies: evolution - galaxies: haloes - 
galaxies: stellar content - (cosmology:) dark matter - cosmology: observations. 



1 INTRODUCTION 

Understanding the origin and evolution of galaxies, in par- 
ticular the most massive, is one of the major challenges 
in modern astrophysics. Many massive galaxies today are 
giant early-type systems; hence the formation of spheroids 
should proceed to a certain extent in locked step with the 
mass assembly. Numerical sim ulations based on th e A Cold 
Dark Matter (CDM) theory (|White fc Reesi Il978| j predict 
the hierarchical mass assembly, i.e., large CDM halos (or 
dark halos) are formed through successive mergers of smaller 
building-block systems. Within formed dark halos, baryons 
dissipate their energy and gravitationally collapse to form lu- 
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minous galaxies. The simplest consequence of the above sce- 
nario is that galaxies are also formed in a hierarchical way, 
which means that the most massive galaxies emerge in the 
last phase of the formation history. Thus the number-density 
measurements of massive galaxies, or luminous red galaxies 
as a proxy of them, in the distant universe are one of the key 
observations in recent studies of cosmology and galaxy evo- 
lution. While there has been growing evidence tha t luminous 
red galaxies are largely i n place at z ~ 1 (e.g., iBell et al.l 
12004; iBorch et al.l [2003: ICimatt i. Daddi. fc Renzinil 1200^ 
IWillmer et al.ll2006l: iBrown et all 120071 : iFaber et al.ll2007| j. 
recently IMatsuoka fc Kawaral ( 2010h found an evidence for 



the hierarchical formation of massive (M* > 10" M ) galax- 
ies since z — 1 based on an analysis of a large galaxy sample 
extracted from the near-infrared images. 

Clustering measurements provide another important 
clue as to the formation and evolution of massive galaxies. 
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In the structure formation paradigm of the A CDM theory, 
dark halos are predicted to cluster together with the cluster- 
ing strengths depending on their masses, in such a way that 
more massive dark halos have larger clustering amplitudes. 
Thus measurements of galaxy clustering can be used to infer 
the underlying CDM distribution, and eventually reveal the 
evolutionary path of galaxies throughout the cosmic mass 
assembly and structure formation history. There have been a 
number of cluster ing measurements of luminous red galaxies 
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20081 ) universe. They found that more lu- 
minous galaxies are generally more clustered, as expected if 
more massive (thus more luminous in red wave bands) galax- 
ies reside in more massive dark halos. Clustering measure- 
ments of stellar-mass selected galaxies are much more chal- 
lenging due to the difficulty in estimating stellar masses of 
a large number of ga laxies over a wide field of sky. Recently 
iFoucaud et al.l l|2010l ) present such a measurement for galax- 
ies with the stellar masses Ah > 10 10 M o at 0.4 < z < 2.0, 
taken from the Palomar Observatory Wide-field Infrared 
Survey, and report the evolution of the stellar-mass to total- 
mass ratios inferred from the clustering strengths of galaxies. 

In this work we present the clustering measurements of 
massive ( M+ > 10 1 Mp ) ) galaxies at 0.2 < z < 1.0, a sample 
built up in iMatsuoka fc Kawaral ((2010), and show a simple 
interpretation of the measured angular correlation function 
by using halo occupation distribution (HOD) models. This 
paper is organized as follows. We introduce our galaxy sam- 
ple and quantify their angular correlation in Section 2. The 
construction of HOD models and the comparison with the 
observations are described in Section 3. The inferred proper- 
ties of underlying dark halos are discussed in Section 4, then 
a summary follows in Section 5. Throughout this paper, we 
adopt the cosmological parameters of Ho ~ lOO/i = 70 km 
s _1 Mpc -1 , n M = 0.3, n h = 0.04, Q A = 0.7, cr 8 = 0.8, and 
n s = 1.0. All magnitudes are expressed in the Vega system. 



2 MEASUREMENTS 
2.1 Sample 

The massive galaxy sample used in this work is taken from 
IMatsuoka fc Kawaral (|2010l ). Below we give a short descrip- 
tion of the sample, while we refer the reader to the above 
paper for the full details. 

The massive galaxies are drawn from 55.2 deg 2 of the 
UKIRT (United Kingd om Infrared Telescope ) Infrared Deep 
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galaxies in the if-band sources are correctly extracted and 
that stellar contamination is negligible. Photometric red- 
shifts have been derived by comparing the galaxy colours 
in the u, g, r, i, z, Y , J, H, and K bands to the spec- 
tral energy distribution templates created from a part of 
the sample with spectroscopic redshifts. The estimated red- 
shift accuracy is 0"Az/(i+ z ) ~ 0.04. Stellar masses have been 
estim ated by fitting the stellar population synthesis mod- 
els of iBruzual fc Charlotl l|2003h to the observed colours, by 
varying the age, metallicity, typical star-formation duration, 
and dust extinction of stellar p opulations. We adopted the 
initial mass function (IMF) of ISalpeterl (|l955l ). The fairly 
robust estimates of stellar mass were obtain ed thanks to the 
inclu sion of near-infrared photometry (e.g., Matsuok a et all 
2008), with the estimated uncertainty o"Aio g M* ~ 0.2. The 
rest-frame U — V colours were calculated by fc-correcting the 
observed r, i, or z-band magnitudes. 

We define the massive galaxies in this work with the 
stellar masses 1O 11O M < AL < 10 12 M Q . The galaxies 
are grouped into four redshift bins, z = 0.2 - 0.4, 0.4 - 0.6, 
0.6 - 0.8, and 0.8 - 1.0. In addition, we classify the galaxies in 
the following two ways. On the one hand, they are divided 
based on their stellar masses into the 10 110 " 11I3 M Q and 
the lO 11 - 5 - 12 ' *^ galaxies. On the other hand, we divide 
them based on their rest-frame U — V colours, into the blue 
(U — V < 1.0) and the red (U— V > 1.0) galaxy populations. 
The total numbers of the subsamples are summarized in 
Table [U 

While the detection completeness is nearly 100 % for 
our sample at z < 1, they could have a significant fraction 
of the contamin ation from less m assive galaxies due to the 
Eddington bias l|Eddingtonll 19131 ). since they are located at 
the steep high-end portion of galaxy stellar mass function. It 



is especially significant for the most massive, 

1O U.5-12.O M0 

galaxies, for which a contaminated fraction could be up to 50 
% based on the Monte-Carlo simulation. Such a contamina- 
tion would generally reduce measured amplitudes of galaxy 
clustering. 



Sky Survey (UKIDSS; IHwrence et all 



Digital Sky Survey (SPSS; lYork et al 



20071) and the Sloan 2.2 Correlation Function 



2000) II Supernova 



Survey (|Frieman et al-llioOSl ) , which enable by far the largest 
survey of massive galaxies to date with robust mass esti- 
mate, reaching back to a time when the universe had only 
half its present size. The target field is a strip region from 
the right ascension l h 15 m to 3 h 6 m along the celestial equa- 
tor, i.e., the declination within ±1.25°. The galaxies have 
been selected from the A"-band images of the UKIDSS Large 
Area Survey Data Release 3 (Warren et al, in prep.), and 
separated from Galactic stars based on their r — z and 
z — K colours. We estimate that more than 95 % of the 



A common statistic characterizing clustering of galaxies is 
the angular correlation function (ACF), w(9). It is defined 
by the probability 5P of finding two galaxies separated by 
an angular distance 9, with respect to that expected for a 
random distribution, i.e., 



8P (x S 2 al [l + w(9)]Sn 1 5Q 2 



(1) 



where <5fii and 8Q2 are elements of solid angle and S ga i is 
the mean surface density of galaxies. The ACF w(6) becomes 
zero for a random or homogeneous distribution. 
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We use the lLandv fc Szalavl l|l993l ) estimator for eval- 
uating the observed ACFs as follows. On a field with Ng 
detected galaxies, the numbers of galaxy pairs Ngg(9) are 
counted as a function of the separation angle 9 ± AO (deg). 
Here we adopt the logarithmic step size A (log 9) = 0.2. On 
the other hand, JVr points are placed at random over the 
observed fields as covered by the real data (galaxies), and 
the numbers of random pairs N R k{9) are counted. In addi- 
tion, the numbers of data (galaxy)-random pairs, Nq-r{9), 
are also counted with the above data and random points. 
Then the lLandv fc Szalavl (|l993T l estimator gives 



w*(9) = 



DD{9) - 2DR(9) + RR{9) 
RR{9) ' 



(2) 



where DD, DR, and RR are the data-data, data-random, 
and random-random pair counts normalized by the to- 
tal numbers of pairs, i.e., DD = Ngg/[^Ng(Ng — 1)], 
DR = N G r/[N g Nb], and RR = N RR /[±N R (N R - 1)]. We 
repeat the calculation 10 times with varying the random 
components. 

The measured w*(9) are known to be biased to lower 
values with respect to the real w(9) due to the boundary ef- 
fect of the observed ar ea. This bias offset is kn own as the in- 
tegral constraint, wic |Croth fc Peebleslll977r ). We estimate 
its amplitude by approximating the ACF with the power-law 
form 



w{9) = w*(6) + wic = A* 



(3) 



In the case of a power-law ACF, the integral constraint is 
given from the random-random pair counts by Wi c/A w = 
E s Nnn(0)9~ s / £ 9 JV RR (0) fe.g.. lRoche et all 19991 ). We as- 
sume the standard value of the ACF slope 8 = 0.8, which 
is found in a wide range of optical an d infrared observa- 
tions of nearby and d istant galaxies fe.g. jBaugh et aL l l 19961 ; 
iRoche fc Ealesl 1 19991 ). While the steeper slopes (<$ ~ 1) 
have been rep orted for nearby early- type galaxies (e.g., 
lLovedav et alJ 1995; Gu zzo et alj Il997f ) and luminous red 
galaxies fe.g.. lBrown et al.ll2008r ). we confirmed that it is suf- 
ficient to adopt a single value of 5 = 0.8 for our sample con- 
sidering the relatively large errors in the ACF measurements 
especially for the most massive galaxies. The estimated am- 
plitudes of the integral constraint are wic ~ 0.04 w(l'), 
which are negligible at most of the angular scales studied in 
this work. 

We divide our strip-shaped field into nine sub-fields 
along the right ascension, each covering approximately 3.0 x 
2.5 deg 2 , and measure the ACF in each of the sub-fields. 
The measurements in the sub-fields are then averaged to 
give the mean estimates of w(9) with the associated errors 
calculated from the field-to-field scatters by the bootstrap 
resampling. The results are shown in Figs. [1] and [2] for the 
different stellar-mass, color, and redshift classes. Note that 
one should take the covariance matrices into account in rig- 
orous analyses of the ACFs, since the measured ACF am- 
plitudes at different angular scales are correlated with each 
other. However, it would have little impact on the conclu- 
sions presented in this work, given that the uncertainties of 
the deduced halo characteristics come mostly from the in- 
adequacy of the current HOD models rather than the mea- 
surement errors (see below). 



3 HALO OCCUPATION DISTRIBUTION 
MODELS 

We investigate the evolutionary link of the observed galax- 
ies to underlying dark (CDM ) halos by using halo occu- 
pation distribution (HOD; see lCoorav fc Shethl l2002. for a 
review) models. HOD modeling is a powerful approach to 
explore galaxy distribution within dark halos, whose for- 
mation and evolution can be predicted by simulations and 
analytic methods. Observations can provide the useful in- 
sights into underlying dark halos by constraining HOD pa- 
rameters in such a way that the models appropriately re- 
produce measured clustering prop erties of galaxies (e.g., 
Zehavi et aill2005l; iLee et all 120061 1. In this work we follow 
Blake. Collister. fc Lahavl (1200 8) for constructing our HOD 
models ( see also I Tinker et al.l l2005). We give a short descrip- 
tion of the models below, while the reader is referred to the 
above papers and the references therein for further details. 



3.1 Model Construction 

In the HOD models we use the h alo mass function n(M) 
presented by Ijenkins et al. ll200ll). We assum e the scale- 
dependent bias following Tinker et al.l (|2005l). with the 
bias f unction b(M) taken from the ISheth. Mo. fc Tormenl 
ll200ll) mode l inclu ding the revised parameters given by 
iTinker et al.1 ((2005). The non-linear dark-matter power 
spectrum is const ructed using the fitting formulae of 



Smith et all (120031). Th e dar k halo density profile of 



Navarro. Frenk. fc White] (|l997t) is adopted, with the char- 



acterizing concentra tion parameter c(M ) calibrated by nu- 
merical simulations (|Bullock et al.ll200ll ). 

The statistical number of a central galaxy within a dark 
halo of the mass M is parametrized as 



N cen (M) = 0.5 



1 + erf 



log (M/M cu t 

Ocut 



(4) 



where erf (a;) is the error function (|Zheng et al.ll2005l ). The 
halos above a threshold mass of M cu t have a single cen- 
tral galaxy and those below M cut have no central galax- 
ies on average, while the transition between the two cases 
is soften by the parameter cr cu t. We adopt CT cut = 0.3, 
which i s the mean value found for lum i nous red galaxies a t 
z < 1 l|Blake. Collister. fc Lahwll200Sl ; iBrown et al.ll2008l ). 
The number of satellite galaxies is described as the power- 
law form 



(5) 



The total number of galaxies populating a dark halo is given 
by 



N(M) = N cen (M)[l + JV Mt (M)], 



(6) 



thus no galaxies are statistically found in the dark halos 
below the threshold mass M cut . We show an example of our 
HOD model along with the adopted halo mass function in 
Fig. El 

A spatial clustering of galaxies £(r) can be divided into 
the two terms, namely, a one-halo term (£ih(f )) arising from 
pairs of galaxies in a same dark halo, and a two-halo term 
(£211 ( r )) arising from pairs of galaxies that reside in different 
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Figure 1. The measured angular correlation functions (diamonds) for the 1O 11O_11 ' 5 M0 (left panels) and the lo 11 - 5-12 Mq (right 
panels) galaxies at z = 0.2 - 0.4, 0.4 — 0.6, 0.6 — 0.8, and 0.8 — 1.0, from top to bottom. The error bars represent the la uncertainties 
of the measurements. The green and red curves show the best-fit HOD models of the case (i) and (ii), respectively (see text); the dotted 
and dashed lines represent the one-halo and two-halo terms, respectively, and the solid lines represent their sums. The thin vertical lines 
mark the comoving scales of 0.5 and 5 h~ x Mpc at each redshift. 
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Figure 2. Same as Fig.[T] but for the blue (left panels) and the red (right panels) galaxy populations. 
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Figure 3. An example of the HOD model N(M) and the halo 
mass function n(M). The HOD model is calculated with the 
parameter values M cut = 10 13 M Q , Mo = 1O 14 M , and /3 = 
1.5. The dotted lines represent the numbers of central galaxies 
N ccn (M) and satellite galaxies N sat (M), and the solid lines rep- 
resent the total number N(M) = N ccn (M)[l + 7V sa t(M)]. The 
dashed lines show the evolving halo mass function at z = 0.3, 0.5, 
0.7, and 0.9, from top to bottom. An actual number of galaxies 
hosted in a halo with the mass M in a given volume scales as 
n(M) x N(M). 



halos: 



eW = [l + Cih(r)]+6h(r). 



(7) 



The one-halo term dominates a correlation function at small 
scales < a few ft -1 Mpc, i.e., less than a typical virial diam- 
eter of large halos, while a two-halo term becomes dominant 
at larger scales. The halo exclusion (which implements the 
condition that halos cannot reside within each other) is con- 
sidered by introducing the truncation mass, above which 
the halos would overlap with each other at the separation 
r, into the calculation of £211 (f)- The truncation mass is 
derived using the ' n'„-matched' approximation described in 



iTinker et al.ll|2005t ). The spatial correlation £(r) is projected 
to the angular correlation function w(8) via 

p(z) 



w(6) = 2 dx 



dx(z) /dz 



du t{\/u 2 + x(z) 2 9 2 ), 



(8) 

where x(z) is the comoving radial distance to the red- 
shift z. The redshift distribution of galaxies, p(z), is es- 
timated by giving the random gaussian errors with the 
standard deviation £>"az/(i+z) ~ 0.04, which corresponds 
to our photometric-redshift uncertaint y, to the interpolated 



numb er-density evolution reported in IMatsuoka & Kawara 



2010) (including the local measurements by I Cole et al 
2001)) multiplied by the comoving volume element, and 
then extracting the true redshift distribution of the sources 
found in the disturbed-redshift bin under consideration. 
However, we find that the resultant HOD parameters pre- 
sented below are not very sensitive to the details of the p(z) 
function. We show the examples of the calculated ACFs with 
the several sets of the HOD parameters in Fig. [3] 

From the above parametrization of the HOD N(M) 
along with the halo mass function n(M) and the bias func- 
tion b(M), we can calculate some useful quantities; the num- 
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Figure 4. The dependence of the predicted ACF on the HOD 
parameters. In each panel the two curves of the same line style 
represent the one-halo and the two-halo terms, the former be- 
ing dominant at small angular scales. The basic parameters are 
(Mcut, cr cut , M , 0) = (1O 13 M , 0.3, 10 14 M Q , 1.5) and one of the 
parameter values is altered in each panel as indicated at the top 
right corner. Note that (T cu t is fixed to 0.3 in the present work. 



ber density of galaxies: 

n ga i = j dM n(M) N(M), 
the effective large-scale bias: 



(9) 



bgai = / dM b{M) n(M) N ^ M \ (10) 



(11) 



Kgal 

the effective mass of the host dark halo: 

NjM) 
n ga i 

and the average fraction of central or satellite galaxies: 
/con =— — f dM n(M) N ccn (M); central, 

"gal J (12) 



Af cff = J dM M n{M) 



/sat =1 - /c Cn ; satellite. 



3.2 Other Models 

Note that the above HOD form is based on a theoreti- 
cal galaxy sample above a certain baryonic mass thresh- 
old, predicted by smoothed particle hydrodynamics (SPH) 
simulations and semi-analytic galaxy formation models 
l|Berlind et all 120031 ; IZheng et alj |2005). While our galaxy 
sample is defined by the stellar-mass bins rather than thresh- 
olds, we use the above models in this work for the follow- 
ing reasons. First, our massive galaxies are located at the 
steep high-end part of the mass function, so that the galax- 
ies above the upper limit of a stellar-mass bin have little 
impact on the observables whether they are included in the 
measurements or not. In fact we find no significant difference 
in the measured number densities and ACFs between the 
M* = 10 110 " 11 - 5 M Q and M* > 10 110 M Q galaxy samples 
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or between the M* = lO 11 - 5 - 12 O M and M* > 10 115 M© 
galaxy samples. Second, we actually find that considering 
an upper cut-off in N cen (M), which is used to reproduce 
the HOD of a bin sample (see below), does not provide the 
better fits to the observed ACFs than our standard model 
does. Third, the adopted model can approximate the HODs 
of various kinds of galaxies with a minimum number of pa- 
rameters (|Zheng et al.ll2005l ), which is essential in exploring 
the halo characteristics of galaxies with relatively large er- 
rors in the ACF measurements, such as our sample. 

We have examined the effects of changing HOD models 
as follows. The standard model is modified by considering 
the variable a cut , fixed f3 and/or the additional cut-off pa- 
rameter Mi of the satellite HOD: 



iV sat (M) 



M — M i \ P 
M 



We have also tried the six-parameter form: 



iVcen(M) = N exp 



log (M/M cen ) 



2a 2 



(No < 1) 



N sat (M) = 



M - Mi 
Mo 



(13) 



(14) 



which is constructed to reproduce the HODs of a sample of 
galaxies in a stellar-mass bin or of young (in terms of stellar 
age) galaxies presented bv lZheng et alj (|2005T ). The essential 
difference of this form from our standard one is that it allows 
an upper cut-off of N CBn (M). In addition, we have replaced 
the whole model construction process (from the choice of 
dark-matter power spectr u m) fro m the above one following 
iBlake. Collister fc Lahavl J2008h to the one following the 
recipe given by Hamana et al.l (|2004l 'l. who adopt the HOD 
form of 



*<«> = {1 M/M ' ) ' 



for M > M min , 
for M < M min . 



(15) 



We hav e also tried a few other halo ma ss func- 
tions JSheth fc Tormenl 1 19991 : iReed et all 120071 : 

iBhattacharva et al. 20ltf) and bias functions 

However, 



i Sheth fc Tormenl 1 19991 : iTinker etafl l20ld) 



we find that none of the above altered models provide 
the better reproduction of the observed ACFs than our 
standard model does. 



3.3 Results 

We constrain the parameters M cut , Mo, and /3 of the stan- 
dard HOD model by comparing the predicted ACFs (w th (8)) 
and number densities (n*^) with those observed (w° bB (0) 
and n°af) through the Markov chain Monte Carlo (MCMC) 
method. The observed nu mber densities are taken from 
iMatsuoka fc Kawaral |201(t ). We evaluate the likelihood of 
the specific model by the x 2 value of the model fit to the 
observed quantities. The ACF fitting is restricted to the an- 
gular scales 8 < 0.3 deg where the effects of the integral 
constraint (i.e., the boundary effect of the observations) are 
negligible. 

We summarize the obtained best-fit HOD parameters in 
Tabled (the former values in the parentheses). In the case of 
the 1O 11I5 " 12 ' O M0 galaxies at 0.6 < z < 0.8, the parameter 
/3 is not well constrained in the plausible range (/3 < 3) and 



we fix it to the best-fit value found for the 10 U-0_n,B M e 
galaxies at the same redshift range (it is justified by the fact 
that we find the consistent /3 values for the two populations 
at other redshifts). The poor constraint on ft is caused by 
the fact that the three HOD parameters correlate with each 
other in reproducing a given ACF, hence the MCMC fitting 
could converge to unrealistic combinations of the parame- 
ter values when a measured ACF is accompanied by large 
uncertainty. In other subsamples, plausible values of /3 are 
obtained; we find ft = 1.5 — 2.0 for all but the blue sub- 
samples, which are in good agreement with those found for 
SPSS luminous red galax i es at 0.4 < z < 0.7 measured by 
IBlake. Collister. fc Lahavl IMPOST ) . 

The ACFs of the best-fit HOD models are overlaid on 
the observed ACFs in Figs. [T] and [2] (green curves). In these 
fittings, the systematic discrepancy between the model and 
measured ACFs are observed at the large angular scales 
where two-halo terms dominate the ACFs. The reason for 
this discrepancy is fairly clear; the observed ACFs are so 
strong that very massive halos are required to reproduce 
them, while the predicted number of such massive halos is 
very small compared to the observed numbers of massive 
galaxies (the further discussion is given below). As a ref- 
erence, we have carried out another set of model fittings 
with the modified requirement that the models should only 
reproduce the observed ACFs, imposing no constraint on 
the predicted number densities of galaxies. We refer to this 
second set of model fittings as case (ii) , and the original fit- 
tings as case (i) hereafter. The results of the former case are 
listed in Table [5] (the latter values in the parentheses) and 
are shown in Figs. [T] and [2] by the red curves. We obtain the 
excellent fits of the ACFs in the case (ii) over the entire an- 
gular scales, while the predicted are smaller than n°^f 
in the most cases. 



4 DISCUSSION 

As a starting point for the discussion, we show the ampli- 
tudes of the measured ACFs at the comoving scales of 0.5 
and 5 ft -1 Mpc in Fig. [5] They are derived by interpolating 
the data points in Figs.[T]and[2]with the best-fit HOD mod- 
els of the case (ii)Q The two comoving scales are chosen to 
represent the clustering within a halo (i.e., one-halo term) 
or between halos (i.e., two-halo term), considering the typ- 
ical diameter of a few /i -1 Mpc for large halos. Apparently 
the ACF amplitudes are larger in more massive or redder 
galaxies at both scales. In other words, we find that more 
mature (hence massive or red) galaxy systems are more clus- 
tered. It is the expected result if the more mature galaxies 
have started stellar-mass assembly earlier within the highly- 
biased region where the structure formation has also started 
earlier. The larger diversity of the ACF amplitudes at 0.5 
h^ 1 Mpc than at 5 ft -1 Mpc reflects the fact that the one- 
halo term is much more dependent on the HOD details of 



1 This does not mean that we accept the halo characteristics 
behind the best-fit HOD models of the case (ii). We just use 
these models as the useful templates to interpolate between the 
measurement points in Figs, [l] and [2] taking advantage of the fact 
that the models provide the excellent fits to the measured ACFs. 
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Table 2. Halo Occupation Distributions of Massive Galaxies. 



Redshift Galaxy 


"gal 
















2 sample 


<log M*> (10~ 4 Mpc -3 ) 


log M cut 


log M 







log M eS 


/cen 


~th / „obs 
"gal / "gal 



0.2 


- 0.4 


10 ii.o- 


" n - 5 M 


11.2 


5.5 ± 0.4 


(12.7, 13.4) 


(13.9, 14.6) 


(1.6, 2.0) 


(1.53, 1.86) 


(13.8, 13.9) 


(0.83, 0.94) 


(0.97, 0.18) 














± (0.1, 0.4) 


± (0.1, 0.3) 


± (0.1, 0.1) 


± (0.01, 0.12) 


± (0.1, 0.1) 


± (0.01, 0.02) 


± (0.06, 0.09) 






-i ^,1 1 r i 

10 11 


12M Mq 


11.7 


0.49 ± 0.08 


(13.7, 13.8) 


(14.7, 15.0) 


(1.6, 1.4) 


(2.13, 2.26) 


(14.0, 14.1) 


(0.92, 0.95) 


(0.98, 0.65) 














± (0.1, 0.4) 


± (0.2, 0.6) 


± (0.2, 0.2) 


± (0.05, 0.18) 


± (0.1, 0.2) 


± (0.01, 0.03) 


± (0.14, 0.39) 






blue 




11.2 


0.9 ± 0.4 


(13.4, 13.7) 


(14.7, 14.9) 


(2.4, 2.4) 


(1.90, 2.10) 


(13.9, 14.0) 


(0.97, 0.98) 


(0.88, 0.48) 














± (0.4, 0.5) 


± (0.4, 1.2) 


± (0.2, 0.4) 


± (0.15, 0.23) 


± (0.2, 0.3) 


± (0.01, 0.01) 


± (0.39, 0.32) 






red 




11.3 


5.0 ± 0.5 


(12.8, 13.3) 


(13.9, 14.5) 


(1.5, 1.9) 


(1.58, 1.83) 


(13.8, 13.9) 


(0.81, 0.92) 


(0.95, 0.25) 














± (0.1, 0.3) 


± (0.1, 0.3) 


± (0.1, 0.1) 


± (0.01, 0.09) 


± (0.1, 0.1) 


± (0.01, 0.02) 


± (0.09, 0.10) 


0.4 


- 0.6 


10 llu 


11 t; i c 
11 -°Mq 


11.2 


3.9 ± 0.3 


(12.9, 13.7) 


(14.0, 15.0) 


(1.6, 1.8) 


(1.73, 2.46) 


(13.7, 14.0) 


(0.86, 0.98) 


(0.93, 0.07) 














± (0.1, 0.2) 


± (0.1, 0.3) 


± (0.1, 0.1) 


± (0.02, 0.10) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.07, 0.02) 






,,,11 K 


1 9 fl n r 

i2 ' U M 


11.6 


0.10 ± 0.02 


(14.0, 13.6) 


(15.2, 14.6) 


(1.6, 2.2) 


(2.97, 2.45) 


(14.2, 14.0) 


(0.97, 0.92) 


(1.03, 4.40) 














± (0.1, 0.6) 


± (0.6, 0.7) 


± (0.2, 0.5) 


± (0.07, 0.29) 


± (0.1, 0.3) 


± (0.02, 0.09) 


± (0.15, 4.69) 






blue 




11.2 


1.1 ± 0.2 


(13.3, 13.7) 


(14.6, 16.3) 


(2.4, 0.9) 


(1.98, 2.33) 


(13.8, 13.9) 


(0.98, 0.99) 


(0.89, 0.32) 














± (0.1, 0.2) 


± (0.1, 0.8) 


± (0.3, 0.4) 


± (0.05, 0.12) 


± (0.1, 0.2) 


± (0.01, 0.01) 


± (0.14, 0.10) 






red 




11.2 


2.9 ± 0.3 


(13.0, 13.8) 


(14.0, 15.0) 


(1.6, 2.0) 


(1.86, 2.61) 


(13.8, 14.1) 


(0.83, 0.97) 


(0.93, 0.07) 














± (0.1, 0.2) 


± (0.1, 0.5) 


± (0.1, 0.3) 


± (0.02, 0.12) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.08, 0.02) 


0.6 


- 0.8 


10 110 " 


" 11 ' 5 M 


11.2 


3.4 ± 0.3 


(12.9, 13.9) 


(13.9, 15.2) 


(1.7, 1.7) 


(1.96, 3.04) 


(13.6, 14.0) 


(0.85, 0.99) 


(0.91, 0.04) 














± (0.1, 0.1) 


± (0.1, 0.3) 


± (0.1, 0.2) 


± (0.01, 0.08) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.07, 0.01) 






10 lls " 


" 12 -°M 


11.6 


0.09 ± 0.02 


(14.0, 14.3) 


(15.0, 15.5) 


(1.7, 1.6) 


(3.34, 3.90) 


(14.2, 14.3) 


(0.96, 0.99) 


(0.88, 0.27) 














± (0.1, 0.2) 


± (0.2, 0.5) 


± (0.0, 0.0) 


± (0.09, 0.23) 


± (0.1, 0.2) 


± (0.01, 0.01) 


± (0.18, 0.13) 






blue 




11.2 


1.0 ± 0.1 


(13.3, 13.7) 


(14.4, 15.6) 


(2.2, 1.1) 


(2.24, 2.70) 


(13.7, 13.9) 


(0.95, 0.99) 


(0.89, 0.25) 














± (0.1, 0.2) 


± (0.1, 0.6) 


± (0.2, 0.4) 


± (0.04, 0.10) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.10, 0.06) 






red 




11.2 


2.5 ± 0.2 


(13.0, 14.0) 


(14.0, 15.3) 


(1.8, 1.9) 


(2.10, 3.32) 


(13.7, 14.1) 


(0.85, 0.99) 


(0.90, 0.03) 














± (0.1, 0.1) 


± (0.1, 0.3) 


± (0.1, 0.2) 


± (0.02, 0.10) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.08, 0.01) 


0.8 


- 1.0 


10 n.o- 


■ 1L5 M 


11.2 


1.9 ± 0.2 


(13.1, 14.0) 


(14.0, 15.4) 


(2.0, 1.7) 


(2.31, 3.71) 


(13.6, 14.1) 


(0.89, 0.99) 


(0.92, 0.03) 














± (0.1, 0.1) 


± (0.1, 0.4) 


± (0.1, 0.2) 


± (0.02, 0.12) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.08, 0.01) 






10 11.5- 


- 12 -°M 


11.6 


0.06 ± 0.01 


(14.0, 14.6) 


(15.0, 16.4) 


(2.1, 1.6) 


(3.72, 5.34) 


(14.1, 14.5) 


(0.98, 1.00) 


(0.93, 0.05) 














± (0.1, 0.3) 


± (0.5, 0.7) 


± (0.3, 0.3) 


± (0.10, 0.45) 


± (0.1, 0.2) 


± (0.01, 0.01) 


± (0.19, 0.04) 






blue 




11.2 


0.59 ± 0.06 


(13.4, 13.9) 


(14.6, 16.9) 


(1.8, 0.8) 


(2.58, 3.42) 


(13.7, 14.0) 


(0.97, 1.00) 


(0.94, 0.15) 














± (0.1, 0.2) 


± (0.6, 0.9) 


± (0.3, 0.5) 


± (0.05, 0.18) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.10, 0.05) 






red 




11.2 


1.3 ± 0.1 


(13.2, 14.1) 


(14.1, 15.2) 


(2.0, 2.3) 


(2.46, 3.87) 


(13.7, 14.1) 


(0.89, 0.99) 


(0.89, 0.03) 














± (0.1, 0.1) 


± (0.1, 0.5) 


± (0.1, 0.3) 


± (0.02, 0.13) 


± (0.1, 0.1) 


± (0.01, 0.01) 


± (0.06, 0.01) 



S 

S3 



Note ( x ) — The two values in the parentheses represent the best-fit HOD parameters in the case (i) (former) and (ii) (latter). 

Note ( 2 ) — All masses are given in units of Mq. 
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Figure 5. The ACF amplitudes at the comoving scales of 0.5 h -1 Mpc (left) and 5 h _1 Mpc (right) measured for the 10 11 ' 0- 11,5 Mq 
(green), the lo 11 - 5-120 Mq (orange), the blue (blue), and the red (red) galaxies. The data points are slightly shifted along the redshift 
axis relative to each other for visibility. 



Table 3. Two cases of the HOD model fitting. 



Case 


Fitting results 




Bias or Halo mass 


0) 
(ii) 


w tb (8) < w° hs (8), n% 
w th (8) = w ohs (8), n*^ 


= n ob f 
gal 

<r n° hB 

< "gal 


lower limit 
upper limit 



Note — See text for the details. 



each galaxy population than the two-halo term is. The two- 
halo term is mainly determined by the spatial distribution of 
host halos rather than those of galaxies within a halo, hence 
its dependence on the HOD parameters are much weaker 
than that of the one-halo term. This is clearly seen in Fig. 
II 

In the HOD model fitting, we consider the two cases 
with the different requirements for the models. They are de- 
fined in H3.3I and are summarized in Table As already 
noted, the observed ACFs are so strong that very massive 
halos are required to reproduce them, while such massive 
halos are extremely rare and could not host the observed 
number of massive galaxies. While a small number of halos 
could host a large number of massive galaxies if the major- 
ity of the galaxies are satellites, such a HOD results in very 
strong ACF amplitudes at small scales, which conflicts with 
the present measurements (hence is not selected in the above 
fittings). The case (i) models provide a comparable number 
of galaxies to the observed by including the host halos whose 
masses are presumably lower than in reality, hence the pre- 
dicted ACFs are weaker than the observed. Therefore the 
derived halo masses and biases in the case (i) should be re- 
garded as the lower limits. The case (ii) takes the opposite 
approach by selecting only the most massive halos in order 
to reproduce the observed strong ACFs, while setting aside 
the majority of the observed galaxies presumably hosted by 
lower-mass halos. Thus the derived halo masses and biases 
should be the strict upper limits. We will discuss this is- 
sue further below. The actual HOD parameter values of the 
galaxies should be somewhere in between those of the above 
two extreme cases, and we conservatively adopt the above 
upper and lower limits in this work. 



In Fig. [6]we show the derived bias factors and halo mass 
estimates as a function of redshift. It is clearly seen that the 
most massive galaxies mark the most biased density struc- 
ture of dark matter distributions, as expected. The threshold 
mass for hosting a central galaxy, M cut , is significantly larger 
for more massive galaxies, which is consistent with the gen- 
eral agreement that the capability of hosting massive galax- 
ies depends strongly on halo mass. The typical mass for host- 
ing a satellite galaxy, Mo, is extremely large (~ 10 14_16 Mq), 
which means that only the most massive halos could host a 
galaxy with the stellar mass exceeding 10 11 Mq as a satel- 
lite. Actually we find that the fractions of central galaxies 
/ccn in the current sample are larger than 90 % in most 
cases. These results are the consequence of the relative weak- 
ness of the observed ACFs at small scales (one-halo terms, 
arising from central - satellite pairs) compared to those at 
large scales (two-halo terms, mainly arising from central - 
central pairs in the current case). While the effective halo 
mass M e ff shows a similar trend to M cut with regard to the 
galaxy stellar mass, the difference of the M e g between the 

10 11.0-11.5 Mq and the 1() 11.5-12.0 Mq gala;xiea are less than 

those of the M cu t. It is because of the negative and steepen- 
ing slope of the halo mass function toward high-mass end; 
increasing M cu t is accompanied by the smaller increment of 
M e ff since there are less and less halos as the halo mass in- 
creases. The stellar-mass to halo-mass ratios are found to be 
M*/M c ff ~ 0.003. We observe no significant difference of the 
ratios between the different stellar-mass classes, while the 
ratios might be very marginally larger for more massive sys- 
tems, which indicates that the fractional mass growth rate 
(M/M) of halos are comparable to or marginally lower than 
those of stellar systems above M* ~ 1O 11 M (i.e., M/M < 
M± /M*). The major mass accretion would have already ter- 
minated for these very massive halos, while the stellar-mass 
assembly of the massive galaxies would also be inefficient 
since they are located near the high mass end of the galaxy 
mass function. 

We note that the derived effective masses of ha- 
los, M cff ~ 10 14 M o , are compar able to the halo 
masses found for galaxy clu sters (e.g.. iLiri fc Mohrl |2004| ; 
iLin, Mohr. fc Stanford! [2004). Considering the high frac- 
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Figure 6. Halo occupation distribution of massive galaxies as a function of redshift (top left: b ga i, the bias factor, top right: M cu t, the 
threshold mass for hosting a central galaxy, bottom left: Mo, the typical mass for hosting a satellite galaxy, bottom right: M e g, the 
effective mass of host halos). The green, orange, blue, and red bars represent the 10 110— 11,5 Mq, the 10 11,5— 12,0 -Mq, the blue, and the 
red galaxies, respectively. 



tion of central galaxies (/ CC n > 0.9), the massive galax- 
ies studied here are presumably equivalent to central 
galaxies of galaxy clusters. Actually, the stellar masses 
of M* > 10 11 M Q are typical of the central galaxies of 
galaxy groups or clusters with the halo masses ~ 1O 14 M 
l|Yang. Mo. fc van den Boschll2008h . which is comparable to 
the current estimates of the effective halo mass. 

Fig. [6] shows that the bias factor clearly decreases with 
decreasing redshift. It implies rapid evolution of dark mat- 
ter distributions while the clustering structures traced by 
the massive galaxies evolve relatively slowly. The halo mass 
M cu t also shows a marginal decline toward the local universe, 
which might suggest that the M* > 10 11 Mq galaxies are 
formed in progressively less massive halos with cosmic time. 
Then it might imply that the fractional mass growth rate of 
stellar systems exceeds those of halos {M/M < M*/M+) for 
the relatively low-mass (M* < 10 Mq) galaxies, a part of 
which eventually evolve into the massive galaxies. 

We observe no clear difference of the HODs between 
the blue and the red galaxies in Fig. [6] It contrasts with the 
measured ACF amplitudes plotted in Fig. [5] where the red 



galaxies clearly show higher degree of clustering than the 
blue galaxies. The difference between the two populations 
disappears in the derived HODs due to the relative scarcity 
of the massive blue galaxies; since the current HOD mod- 
els predict that both populations are mostly central galaxies 
and assume that the probability of hosting a central galaxy 
is a monotonically-increasing function of halo mass, the nu- 
merous red central galaxies require relatively low-mass halos 
as their hosts compared to the scarce blue galaxies. Thus the 
derived halo masses (and biases) of the red galaxies become 
relatively low, which cancels out the higher ACF amplitudes 
measured for the red galaxies (indeed the larger halo masses 
or biases are predicted for the redder galaxies in the case 
(ii) which imposes no constraints on the predicted galaxy 
number density). 

However, the HOD assumption that the probability 
of hosting galaxies is a monotonically-increasing function 
of halo mass could be too simplistic particularly for mas- 
sive blue galaxies. It is possible that galaxies in mas- 
sive, grown-up halos have exhausted most of their gas 
and have systematically lower level of star formation. 
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In addition, recent studies suggest star-formation quench- 
ing mechanisms of massive galaxies as a feedback pro- 
cess of galaxy mass ass embly, such as the onset of ac- 
tive g alactic nuclei (e.g., Silk fc R oes 1998; Granato et al.l 



12004 ISpringel, Pi Matteo. fc Hernquistl bOOa T These pro- 
cesses could systematically reduce the number of massive 
blue galaxies in massive halos, in a way dependent on halo 
characteristics and galaxy star-formation and its quenching 
physics. We also note that fitting the ACFs of the blue and 
red galaxies separately with the current HOD models might 
be problematic, since the halos hosting a red central and 
blue satellites, or a blue central and red satellites, cannot be 
properly reproduced. These situations might be mimicked 
by the HOD forms such as those given by the equation (|14[) . 
which allows an upper cut-off in JV cen and reproduces th e 
HODs of 'y our rg' galaxies presented bv lZheng et al.l |2005), 
but we do not obtain improved fits of the ACFs with this 
form as compared to those with our standard model. How- 
ever, since the HODs of the massive blue galaxies could be 
significantly altered by considering the more numerous red 
populations, the central fractions of massive blue galaxies 
deduced above are rather tentative. Further investigation 
of this issue is beyond the scope of this paper, and would 
require much more precise clustering measurements of mas- 
sive blue galaxies. If the observed scarcity of the massive 
blue galaxies is due to the star-formation quenching physics 
rather than their need for very massive host halos, then one 
might be able to conclude, simply from Fig. [5] that massive 
blue galaxies reside in less massive halos than the hosts of 
massive red galaxies. Again, it would be the natural result 
if redder galaxies have started stellar-mass assembly earlier 
within the region where the halo mass assembly has also 
started earlier. 

Finally, we comment on the discrepancy between the 
HOD model predictions and the observed quantities. As 
stated above, the observed ACFs are so strong that very 
massive halos are required to reproduce them, while such 
massive halos are extremely rare and could not host the ob- 
served number of massive galaxies. It is unlikely that most of 
the massive galaxies are hosted by a small number of such 
massive halos as satellites, since it would result in much 
stronger ACFs than observed at small angular scales. The 
discrepancy between the observed and model ACFs in the 
case (i) tends to increase with increasing redshift (see Figs. 
[T]and[2|. In this regard, we note that a similar discrepancy 
is observed in the analysis of distant red galax ies at z ~ 2.3 
by_ Tinker. Wechsler. fc Zhenel (|2010h (see also lQuadri et al.l 
(2008)). The systematic difference between their model and 
observed ACFs at large angular scales is roughly a factor 
of two (while within the 1 a uncertainty of their measure- 
ments), whic h is comparab l e to t hose found in our sample. 
As raised by lOuadri et al.l (|2008l ). the underlying cause of 
the problem might be that the current HOD models are too 
simplistic. While the current models assume that a num- 
ber of galaxies hosted in a halo depends solely on the halo 
mass, other halo characteristics such as the mass accretion 
rate could affect the galaxy formation and alter the observed 
properties of galaxies within the halo. We also note that fit- 
ting the clusterings of different galaxy populations (e.g., blue 
and red galaxies) at the same time may improve the situa- 
tion, while such an analysis would significantly increase the 
number of HOD parameters and require much precise data 



for clustering measurements. Furthermore, we can raise the 
possibility that the evolution of halo mass function and/or 
bias function is not well understood; our results might point 
to the earlier emergence of massive halos (i.e., larger num- 
bers of massive halos at high redshifts) than predicted in 
the current halo models. 

It is also worth pointing out that the observed ACFs 
are apparently larger than the fitted HOD models at the 
smallest angular scales (6 ~ 0.001 deg) . These close pair s 
may represent interacting galaxies (e.g.. lRoche et al]|l999l ). 
Reproducing such a feature by the HOD models requires 
the detailed description of galaxy interaction within a halo. 
However, our focus here is the global characteristic of the 
host halos and such an investigation is a subject of future 
papers. 



5 SUMMARY 

In this work we present a clustering analysis of ~60,000 
massive (stellar mass M±> 10 Mq) galaxies taken from 
iMatsuoka fc Kawaral (|2010l ). The galaxies have been ex- 
tracted from 55.2 deg 2 of the UKIRT Infrared Deep Sky Sur- 
vey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II 
Supernova Survey, and are nearly complete to z — 1 with ro- 
bust estimates of photometric redshift and stellar mass. We 
classify them based on stellar masses 

or 1O U - 6 " 12 '°M ) and rest-frame colors (blue: U-V < 1.0, 
red: U — V > 1.0), in order to reveal the difference in the 
spatial distributions of different galaxy populations. 

The angular correlatio n functions (ACF s ) of th e galax- 
ies are quantified using the lLandv fc Szalavl ljl993h estima- 
tor, and we find strong clustering detected for all the sub- 
samples. In order to interpret the measured ACFs, we con- 
struct halo occupation dis t ribu tion (HOD) models following 
iBlake. Collister. fc Lahavl IMPOST) and compare the predicted 
ACFs with the observed. Our major findings are as follows. 

(i) The clustering amplitudes clearly show a systematic 
trend regardless of the measured redshift, in which the 
most massive, 1O 11,5 ~ 12 '°M galaxies have the strongest 
clustering and the blue galaxies have the weakest cluster- 
ing. It would be the natural result if more mature galaxies 
have started stellar-mass assembly earlier within the highly- 
biased region where the halo mass assembly has also started 
earlier. 

(ii) The bias factors and halo masses are systematically 
larger for the galaxies with larger stellar masses, which con- 
firms that the capability of hosting massive galaxies depends 
strongly on halo mass. The stellar-mass to halo-mass ratios 
are M*/ M e g ~ 0.003, with no significant difference observed 
between the different stellar-mass classes. 

(iii) The inferred halo masses of M eS ~ 1O 14 M and the 
high central fractions (/ CC n > 0.9) indicate that the mas- 
sive galaxies in the current sample are equivalent to central 
galaxies of galaxy clusters. 

(iv) The bias factor decreases with decreasing redshift, 
which implies rapid evolution of dark matter distributions 
while the clustering structures traced by the massive galax- 
ies evolve relatively slowly. The halo mass M cut also shows 
a marginal decline toward the local universe, which might 
suggest that massive (M* > lO n M ) galaxies are formed 
in progressively less massive halos with the cosmic time. 
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(v) If the observed scarcity of massive blue galaxies is due 
to the star-formation quenching physics rather than their 
need for very massive host halos, then one might be able to 
conclude that massive blue galaxies reside in less massive 
halos than the hosts of massive red galaxies. 

(vi) The observed ACFs are so strong at large angular 
scales that very massive halos are required to reproduce 
them, while such massive halos are extremely rare and could 
not host the observed number of massive galaxies. It might 
point to the inadequacy of the current HOD models or a lack 
of knowledge about the evolving halo mass function and/or 
bias function. 
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